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Abstract 

We examine the prospects of detecting an analogous process of neutrinoless double 
beta decay at a neutrino factory from a high energy muon storage ring. Limits from 
LEP experiments, neutrinoless double beta decay as well as from global fits have to 
be incorporated and severely restrict the results. We investigate what limits on light 
and heavy effective Majorana neutrino masses can be obtained and compare them 
with existing ones. Discussed are also contributions from right-handed neutrinos 
and purely right-handed interactions. Other "new physics" contributions to the 
same final state might produce large event numbers. 
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1 Introduction 



The physics potential of a muon storage ring is rich and exciting. Especially the option 
of using the neutrinos from the /i decay gathered much attention Typically the main 
focus lies in long baseline oscillation experiments with source-detector distances from 
730 up to 10000 km. This development is driven by the urge to find out about oscillation 
phenomena in more detail and to gain additional information, be it about CP violation, 
the sign of Am^, the size of \Ue3\ or the existence of sterile neutrinos. 
An additional option is the usage of a detector directly at the storage ring site. Neutrino 
interactions of up to 10^^/yr provide the possibility of high precision experiments regarding 
CKM matrix elements, structure functions, electroweak parameters, charm physics or other 
phenomena, see for some possibilities. As in any other new experiment, new physics may 
lurk in the results. In the light of recent developments in oscillation experiments, effects of 
massive neutrinos are hot candidates. Evidence for massive neutrinos and therefore physics 
beyond the Standard Model (SM) comes from the up-down asymmetry of atmospheric 
muon-neutrinos, the deficit of solar neutrinos and the LSND experiment. See |^ for more 
complete surveys. For example, the see-saw mechanism might connect the very light 
known neutrinos to heavy neutrinos, which are usually assumed to be of Majorana nature. 
Majorana particles can show their presence not only by being directly produced, but also 
via indirect effects stemming from their B — L violating mass term. The best known 
example for such a process is neutrinoless double beta decay (Oz//5/5) 0, which results in 
limits on the effective electron neutrino Majorana mass (rriee)- The complete 3x3 matrix 
of (light) effective Majorana masses is defined as 

{friap) = \iU diag(mi,m2,m3)[/'^)„^| 
= \2_^miUaiUf3i\ < 2_^m.i\UaiUpi\ with = e, /x, r, 

where the sum goes over the mass eigenstates m^. Conversely, the "inverse effective mass" 
is defined as 
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1 



^ — Uc^iUpi < — \UaiU(ii\ with a,/? = e, /x, r. 



(2) 



The sum over i is not the same in Eqs. and (0): For {map) it goes over all "light" 
mass eigenstates and in {j^—^) over all "heavy" states. The attribute "light" or "heavy" 
depends on the energy scale of the process one considers to obtain information about the 
respective element. Note that we indicated that the sum goes up to a number greater than 
3, it is however also possible that only one additional very heavy neutrino exists. Apart 
from theoretical prejudices, a priori we do not know how many there are. 
The latter matrix might seem somewhat artificial; its form comes from the fact that heavy 
Majorana neutrinos force cross sections or branching ratios (typically processes in analogy 
to 0uf3j3) into a mass~^ behavior. The knowledge of the elements of the matrices is rather 
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poor, of course with the exception of (mee) and (jj^)- At a neutrino factory, the following 
processes (see Fig. |I]) can be used to gain information about the other elements: 

^iy\ N where / = e, /i and a, /3 = e, /i, r. (3) 

Here X denotes the hadronic final state. Because of very general arguments §, the pure 
observation of this process guarantees a. B — L violating Majorana mass and in connection 
with SUSY a B — L violating sneutrino mass term. This connection is depicted in Fig. ^ for 
the non-SUSY case. The precise determination of the mass of the intermediate Majorana 
neutrino will be very difficult, however, even the demonstration of Majorana mass terms 
will be an exciting and important result, since different models predict different mass 
matrices. In some models {rriee) is zero and therefore the only direct information about 
the mass matrix might come from neutrino oscillations. This complicates the situation, 
since only mass squared differences are measured and the additional phases induced by 
the Majorana nature are unobservable. Other experiments or cosmological arguments give 
total mass scales but the precise matrix is highly nontrivial to find . Thus, information 
about entries in {rriais) is very important. Similar arguments hold for the existence of heavy 
Majorana neutrinos. Note that the see-saw formula predicts their mass to lie in the 
range 

mN^^^lO\..WGeV, (4) 

where is a charged lepton or quark mass (i.e. electron to top quark) and the mass 
of a light neutrino (10~^ ... 1 eV as indicated by oscillation experiments). It turns out that 
the highest cross section of process (|^) is obtained for the lower region of this mass range. 
In addition, if a -B — L violating process is detected, it is helpful to know if the "mildly 
extended" (i.e. just additional Majorana neutrinos) SM can provide the signal or another 
theory, such as SUSY has to be considered. 

The paper is organized as follows: In Section ^ we discuss some properties of process (0) 

and its application to neutrino factory kinematics. We review in Section ^ the status of 

direct experimental limits on {rriap)- For the first time we give — using HERA data — 
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bounds on elements of ) other than (i^^) and examine what new limits might be 
accomplished for different neutrino factories. It turns out that for muon energies higher 
than 500 GeV physically meaningful limits on {rriais) can be obtained. Then we summarize 
limits on heavy Majorana neutrinos and their mixing with SM particles. Regarding the 
prospects of detecting events from process (|^) we apply in Section ^ all these limits, which 
severely restricts the results. 

2 The process and a neutrino factory 
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2.1 Kinematics 



The Feynman diagram is shown in Fig. |l|, the calculation is straightforward and described 
in more detail in [l^. For one eigenstate rrii it was found: 



IMWuiq^l-a+a+q') 
1 



1 



w) 



{Pl ■ P2) 



7^ 



{ki ■ k3){k2 ■ ki 



■{{k2 ■ k3){ki ■ ki) - {ki ■ k2){k3 ■ k^) - (ki ■ k3){k2 ■ k^)) 



{Qi-mi){q2^ - mi) 

(5) 

Here q2 denotes the momentum of the Majorana neutrino in the crossed diagram, which has 
a relative sign due to the interchange of two identical fermion lines. In addition one has to 
include a factor | to avoid double counting in the phase space integration. Scattering with 
an antiquark and with an antineutrino is equivalent to the following simple replacements: 



\M\''{uiq^ I- a+a+q') = \M-+\^ = \M—Wp2 ^ k^), 



\M\\Viq^V 



a a q') 



a a q') 



\M+-\ 



\M\ 



(6) 



Of course, the two leptons from the intermediate ^WW — ^ af3" diagram do not have to be 
of the same flavor: An interesting statistical effect P] occurs when one considers the relative 
difference between, say, the /i/i and the /le final state (mass effects play no role for e and /i): 
First, there is no phase space factor ^ for the latter case. Then, there is the possibility that 
an electron is produced at the ("upper") ulW vertex or at the ("lower") qq'W vertex. Both 
diagrams are topologically distinct and thus have to be treated separately. This means, 
four diagrams lead to the /xe final state, whereas only two lead to the /i/i final state. We 
see that there is a relative factor 4 between the two cases. Note though that now the 
interference terms are added to the two squared amplitudes since there is no relative sign 
between the two. This reduces the relative factor to about 3. 



The details of our Monte Carlo program are given in |10], now an integration over the 
incoming neutrino energy spectrum has to be included. A simple phase space calculation 
for the yU — > eueiy^^ decay gives for the normalized distribution in the lab frame: 



dN 



dN 



E\ 



(7) 



The maximal neutrino energy is E^ and the mean value is (E,^) = 7/10 (3/5) i?^ for z/^ 



In Figs. to |5| we show the total cross section for the reaction Uf^N ^ a'^f3'^X with the 
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/i/i, /ir and rr final states for three different muon energies. For the above given statistical 
arguments to be valid, the mass of the final state leptons has to be negligible. This is the 
case for muon energies higher than about 3 TeV, as can be seen from the figures, where 
the /ir final state is the leading signal for = A TeV. We set Uai = 1 in order to show 
the mass dependence of the signal; the slope of the curves is easily understood from the 
two extreme limits of 



^2 C for mj < 



(q -rrii) 1 ^ for > q 



2 



where q is the momentum of the Majorana neutrino. 

For a yU~yU "''-collider or a muon storage ring four different signals are possible (corresponding 
to incoming i/^, !7g, 17^ or z/g). Fig. ^ shows that the muon neutrino from the /i~ decay gives 
the highest cross section. In this figure we plotted the interesting area of the mass range 
given by the see-saw formula (§) and applied also the limits on heavy neutrino mixing as 
explained in Section Finally, we give in Fig. |^ the cross section for u^N 
with two possible other realizations, namely via an intermediate right-handed Majorana 
neutrino and via right-handed interactions with a WO^-mass of 720 GeV, the current lower 



bound [ITJ]. For the latter case, the matrix elements are identical whereas for the former 
one has to make the replacement {p2 ^ ^4, Pi ^ fci) in Eqs. (|^) and d). It can be seen, 
that a left-handed heavy neutrino gives the highest contribution. Of course it is possible 
that all these realizations contribute and thus interfere. 

We checked the dependence of the results on oscillation parameters by integrating over 
two-flavor formulas. Even for E^^ = 50 GeV, a detector-source distance of 1 km and 
LSND-like values of Am^ (~ 0.1 eV^) and sin^ 29 (~ 10~^) the relative suppression of the 
signal was not more than O{10~^). Inserting typical parameters of atmospheric or solar 
experiments has even less effect. 



2.2 Neutrino factories 

Several proposals for a muon storage ring have been discussed, the number of expected 
neutrino interactions differs. The formula used for the luminosity in units of cm~^ s~^ is 



jC = NaJN^I, (9) 

where Na is the Avogadro number, A*"^ the number of muons injected in the ring per 
second, / the fraction of the collider ring occupied by the production straight section and 
I the mass depth of the target in g cm~^. Typical numbers are / = 0.02, / = 1000 g 
cm~^ and A^"^ = 10^^ . . . 10^^ s~^. Let us be optimistic and assume A^ = 10^^ s~^ with a 
"year" of 10^ s running time. With this parameter set one gets C ~ 10'^^ cm~^ s"^. The 
neutrinos from the /i-decay will all end inside the detector since their opening angle is just 
^ ~ 1/7^ = rrin/En. Typical distances between detector and muon ring are 10^ to 10^ m, 
discussed energies go up to 10^ GeV. A complete scope of all possible options is not our 
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aim. If definite plans for experiments and machines are made, our results can easily be 
rescaled with the help of relation (P). 



3 Limits on neutrino parameters 



1 



As expected, Oz//3/3 provides us with the best limit of all entries in {rriap) and {jyi^^ 



Recently, other elements of the mass matrix were investigated and for the first time limits 
on the T sector of {niap) were given |TE[. The process discussed was e'^p Ve a^[3^X at 
HERA and gave bounds on (me,-), (^^t) and {m^T-). In improved limits on {m^^) (via 
— s> TT^fi^fi'^) and {mf.^) (/i~-e^ conversion on titanium) are given. Together with the 
Oi//3/3 limit |T5[ the current situation is as follows: 



{map) 



< 



/ 2 ■ 10-1° [15] 



v 



1.7 (8.2) ■ 10-2 |T4| 
5.0 ■ 10^ P 



4.2 ■ 10^ [0 \ 
4.4 ■ 10^ [|1|] 



GeV. 



(10) 



2.0-104 ^ J 



There is a spread over 14 orders of magnitude. For {nif.^) two values are given, depending 
on the spin configuration of the final state protons. Note that for all entries except for the 
ee element the limits lie in the unphysical region, e.g. for a (mg^) = 4.2 ■ 10^ GeV the cross 
section is proportional to m"^ and not to m? as assumed to get the limit. Improvement 
on most values might be expected from B decays |[T^, . 

For {if^) a limit from Oz//5/5 exists |T^, for which a heavy neutrino has ^ 1 GeV. 



Beside neutrinoless double beta decay there are other ways to get information about heavy 
neutrinos: The LEP machine can produce heavy neutral leptons via e^e~ NN, the 
most stringent limits come from the L3 collaboration |19| ; they exclude masses below 70 to 



80 GeV, depending on the charged lepton they couple to (e, /i or r). On the other hand, 
if heavy neutrinos mix with their light SM counterparts, they should alter the results for 
II decay, u scattering and so on. Global fits then limit the mixing parameters [^, in total 
the limits read: 

Elf/eiP < 6.6-10-^ 



J:\Ur^\^ < 1-^ 



10-3, 
10-2, 



mi > 81.8 GeV, 
rrii > 84.1 GeV, 
m,- > 73.5 GeV. 



(11) 



In Section |4.2| we will discuss how these bounds restrict the possibilities of observing the 
OuPP analogue at a neutrino factory. Before that we apply the procedure from again 
to the HERA data and gain limits on the other elements of ). Here, heavy neutrinos 
must have rrii ^ 100 GeV. The matrix reads: 

/ 1.1 ■ 10-® 5.4 



< 



8.4 



10-3 
10-3 



8.6-10-3 \ 
9.0-10-3 
0.1 



GeV 



-1 



(12) 
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Now there is only a spread of 7 orders of magnitude. All non-ee entries are unpliysical, 
e.g. for the /z/i element one gets with the bound from Eq. ( |ri|) 

\U -P 

rrii > ' f ' ~ 0.7 GeV. (13) 

All limits on the same quantities for a right-handed Majorana neutrino with the usual 
couplings to the SM particles lie in the same order of magnitude. Table |I] shows what 
limits could be achieved for luminosity per year given by the parameter set after Eq. (P). 
The improvement would be tremendous and already for muon energies higher than 500 
GeV the bounds on {niap) lie in the physical region: The limit on {m^^) is about 3 GeV, 
where the slope of the cross section is still rising, i.e. proportional to mf. For {jf^) the 

situation is different: For = 500 GeV the limit on {j^—) is about 0.2 GeV, which 
translates in rrii ^ 0.1 GeV, which is a light neutrino, i.e. in that region is a oc mf. Here, 
energies around 10 TeV are required to get physical meaningful values. 



4 Detection of the process 
4.1 Experimental considerations 

Because of the smallness of the cross section of the process discussed here, one might 
ask if SM processes exist, which fake the signal. A discussion of this kind has already 
been done for trimuon production in uN scattering at previous fixed target experiments, 
both experimentally and theoretically [^. This trimuon production has a ( — h +) 
signature. Due to the principle creation of conventional neutrino beams by using pion- 
and kaon-decays, there is always a "17^-pollution" in the beam which can give a ( — h +) 
signal through muon pair production, be it radiatively or in the hadronic final state via 
e.g. vector meson production. These effects exist on the level of about 10~^ of the total 
observed charged current events. Kinematical cuts to suppress this background, e.g. using 
the invariant mass or angular isolation have been developed. For previous experiments 



however, it was found |jTO[ that for trimuon production via Majorana neutrinos the signal- 
to-background ratio is far too small. However, for a muon storage ring we know exactly 
what neutrino flavor is coming in and thus in the case of /i~-decay there is no SM process 
to give a ( — h+) event. The only exception is z/^ — > fi~e^fi~^X which might be faked by 
a I7e CC event with ^~ production in the jet or via bremsstrahlung. However, as we will 
show below, final states with electrons cannot be expected due to the severe limits from 
Ou/Sp. Because of the z/^z/g or u^u^ structure of the beam, ratios between observed types 
of events (coming from each neutrino species) could be used to establish a signal. Also 
polarization of the muon beam could be useful because it allows to change the neutrino 
spectra and therefore the event ratios in a predictable way. Possible channels involving r 
leptons in the final state might be investigated by topological and kinematical methods as 
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used by CHORUS and NOMAD. Finally, it could even be possible to obtain information 
about CP violation by comparing event numbers from different channels. 
Up to now we ignored in this section effects of neutrino oscillations. An incoming Vg, could 
oscillate into a and create via the aforementioned processes a ( — h +) signal. The 
relevant oscillation parameters are now given by atmospheric (Am^ ~ 10~^ eV^) and 
CHOOZ (sin^ 26 ^ 0.2) data. Integrating the CC cross section of the over a two-fiavor 
formula and taking also into account the factor 10~^, yields numbers smaller than the ratio 
of process (|]) with the (CC + NC) cross section by at least one order of magnitude, even 
for a L = 1 km and = 50 GeV option of the experiment. 

For the other final states there is no SM background: Typical events with additional leptons 
are production of gauge bosons, which however are always accompanied with neutrinos or 
extra jets and thus in principle distinguishable. 



4.2 Is it observable? 

Unfortunately, the bounds on neutrinos and their mixing severely restrict the prospects 
of detecting a signal from B — L violating mass terms at the discussed experiments: For 
example, let us consider a 4 TeV muon source (be it a collider or just a storage ring) and the 
/x/x channel. For the moment, we stick to one rrii. The maximum cross section is achieved for 
a mass eigenstate of about 10 GeV, a^^{mi ~ 10 GeV) ~ 10~^° b. A few years of running 
with C ~ 10'^^ cm~^ s~^ per year could establish an observation. However, for the minimal 
allowed mass of 84.1 GeV, the cross section reduces to a^^inii = 84.1 GeV) ~ 2.0 • 10~^^ 
b, which is then further suppressed by the Uf^i limit to 7.3 ■ 10~^^ b. Roughly the same 
number holds for the ee channel, and for the e/z channel about three times this number. 
However, now the value from Eq. (|T2|) comes into play: Assuming one mass eigenstate 
of 81.8 GeV one gets |?7eip < 9 ■ 10"^ resulting in aeeimi = 81.8 GeV) ~ lO^^s b! 
The cross section stays constant till = 6 ■ 10^ GeV and scales with m^'^ for larger 
masses. For the e/j, and er channel the cross sections are — 2 ■ 10~^^ b/mj[GeV] and 
CTer — 8 • 10~^^ b/mj[GeV], respectively. Thus, the electron final states of process (|) 
provide no real chance for observation. The importance of using the heavy neutrino bound 



from OuP/S was also discussed in p6|, ^| 



Now we investigate possible event numbers: to be independent on the concrete values of 
the experimental parameters we calculate the charged and neutral current cross section by 
integrating over the energy spectrum with the GRV 92 and 98 |^ parton distributions 



including c and b quark contributions. With this we give the maximal ratio (i.e. applying 
all limits of Eq. ([TT|) for the cross section) of the process (^ as shown in Table |^. We 
considered only the muonic and tauonic final states and took for the /ir channel the value 
rrii = 84.1 GeV. The last column gives the number of (CC + NC) events with the 
parameter set given after Eq. (^. Only the highest discussed energy provides a chance 
for observation. However, the realization of this kind of machine remains doubtful, but 
might be realized in a form of a new high-energy physics laboratory In Ref. it 



has been shown that the bounds from Ou(3P on heavy Majorana neutrinos can be evaded if 
one assumes rather baroque models, which seem highly unnatural and are not considered 
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here. 

Though the numbers are no reason to be overoptimistic, the same final state we discussed 
might have contributions from other channels as the ones plotted in Fig. |^ For Oz//5/5 many 
limits on beyond-SM parameters were derived, see for a review. A simple estimation 
shows the power of such a neutrino factory: For a 4 TeV energy and a 100 keV neutrino, 
the cross section is about 10^^^ b. Other contributions might not need a helicity flip and 
are thus larger by roughly a factor of {mj^/E^Y — 10^^, where is the energy of the 
Majorana neutrino. With the mentioned 10^^ cm~^ s~^ luminosity we would have 10^ 
events per year, a "new-physics factory" . 



4.3 Outlook 



Will the situation change with future improved mass and mixing limits? First, the high 
number of neutrino interactions might have impact on the global flts for the mixing matrix 
elements. However, the limits cannot be expected to be improved by factors larger then 
(9(1). The LEP bound on the neutrino mass from |]Ty| correspond to about 40 % of the 
used center of mass energy of 189 GeV, for simplicity we can assume that this will hold 
also for the upgrade energies. As other machines are concerned, at LHC [^, or HERA 
30| investigation of masses of a few 100 GeV might be possible, but mostly only the 
electron channels were considered. Applying the Oz//3/3 limit on heavy neutrinos (which 
has not been done in the analyses) on those results reduces the mass limits in those works. 
For the muon channels the results are signiflcantly lower ^9] or, as for HERA, not yet 
discussed. Regarding NLC, the pair production e^e~ NN produces too small event 
numbers |^l]]. Recently it was shown in |^2|] that for ^/s ^ 500 GeV the "indirect" process 
e^e~ — > ve^W^ might probe Majorana masses up to the center of mass energy. The same 
holds for the e/x option of future colliders via e^yU^ ul^W^ and also Majorana 
neutrino pair production will be observable [^. Processes such as e~7 —>■ u^a' j3~W~^ |^| 
with a, P = fi, T or e~e~ — > [BH] might also evade the OuPP constraint but require 



center of mass energies in the same region as the ones discussed here. 



5 Conclusions 

To conclude, we did a full analysis of the analogue of neutrinoless double beta decay at a 
neutrino factory whilst applying several experimental limits. The use of a detector right at 
the muon storage ring provides a very large number of neutrino interactions and for current 
and future mass limits the signals are perhaps observable at very high muon energies. 
Furthermore, if observed at lower energies, it is important to know how heavy Majorana 
neutrinos with SM couplings contribute to the events. The cross sections are very small 
but at least unaffected by oscillation phenomena. The limits on the effective Majorana 
mass matrix can be pushed down to physical values even from energies of E^j, = 500 GeV 
on. For its "inverse" however, energies higher by a factor of more than 20 are required. 
Signatures of the discussed events might be the only chance to flnd out about Majorana 
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mass terms since most other related B — L violating processes suffer from tiny ratios to 
the respective Standard Model events. Information on future limits on Majorana masses is 
mostly found in works concentrating on electron final states and thus unimportant when 
one incorporates Oz//5/5 bounds on heavy Majorana neutrinos. Anyway, our results should 
not change dramatically even for new LEP limits on direct production or modified global 
fits. Contributions to the same final state without intermediate Majorana neutrinos (e.g. 
SUSY particles) however are a very realistic candidate for observation and the prospects 
for this will be addressed in future works. The process considered in this paper is then the 
relevant background signal. 
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Table 1: Obtainable limits for (m^/j) and {j^—^) (in GeV and GeV-^) for different muon 
energies in GeV. It holds (mge) — (w^^j^) — 1/3 (me^) and {j^) — iff^) — 1/3 ■ 
For the number of events the parameter set given in Section |2]^ is used. For {j^—^) the 
limits from on \Uai\'^ from Eq. [Tl| have been applied. 
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Table 2: Maximal ratio of u^N ^ n^a'^jS^X and sum of CC and NC for for different 
final states and muon energies in GeV. Indirect bounds on mixing matrix elements are 
applied. The last column displays the expected number of (CC + NC) events from 
with the optimistic parameter set given in Section |2.2| . 
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Figure 1: Diagram for v^q — > ^~a^j3^q'. Note that there is a crossed term and ioi a ^ j3 
there are two possibihties for the leptons to be emitted from. The leptonic part also can 
be replaced by the corresponding other neutrino species. 
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Figure 2: Connection between Majorana mass term of v^t and and the existence of 
process v^N — > //~q;+/3+X. 
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Figure 3: Total cross section for Vp^N jj, a^f^^X as a function of the Majorana mass 
for a ii~ energy of 50 GeV. No hmits on Uai are apphed. 
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Figure 4: Total cross section for v^j^N ji q;+/3+X as a function of the Majorana mass 
for a ii~ energy of 500 GeV. No limits on are applied. 
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Figure 5: Total cross section for v^N — > /i a+/?+X as a function of the Majorana mass 
for a energy of 4 TeV. No limits on Uai are applied. 
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Figure 6: Total cross section for UiN IjijiX for v^, u^, and function of the 

Majorana mass for a /i energy of 50 GeV. The hmits on Uai are apphed. 
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Figure 7: Total cross section for u^N — > n'a'^P'^X as a function of the Majorana mass for 
a energy of 50 GeV and different possible reahzations of the process. A?" is a left-handed, 
Nji a right-handed Majorana and Wr denotes the process with a right-handed W boson 
and a left-handed Majorana. No limits on Uai are applied. 
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